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ABSTRACT 

The near-infrared integral field spectrograph SINFONI at the ESO VLT opens a new 
window for the study of central supermassive black holes. With a near-IR spatial 
resolution similar to HST optical and the ability to penetrate dust it provides the 
possibility to explore the low-mass end of the M,-<j relation (a < 120 km s~^) where 
so far very few black hole masses were measured with stellar dynamics. With SINFONI 
we observed the central region of the low-luminosity elliptical galaxy NGC 4486a at 
a spatial resolution of ~ 0.1'^ in the K band. The stellar kinematics was measured 
with a maximum penalised likelihood method considering the region around the CO 
absorption band heads. We determined a black hole mass ofM, = (1.25t^:^^)xl0^ M© 
(90 % C.L.) using the Schwarzschild orbit superposition method including the full 2- 
dimensional spatial information. This mass agrees with the predictions of the M.-cr 
relation, strengthening its validity at the lower a end. 
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1 INTRODUCTION 

Studies of the dynamics of stars and gas in the nuclei of 
nearby galaxies during the last few years have established 
that all galaxies with a massive bulge component contain 



a central supermassive black hole (SMBH; Richstone et al. 
|1998| [Bender k Kormendy |[2003t . Masses of these SMBHs 
(M«) are well correlated with the bulge luminosity or mass, 



respectively, and with the bulge velocity dispersion cr ( Geb 



hardt et al."2000b| [Ferrarese k Merritt||20QQl [Kormendy 



Richstone 1995). There are, however, still a lot of open ques- 



tions in conjunction with the M,-cr correlation, among them 
the exact slope, its universality and the underlying physics. 
A more precise knowledge of the behaviour of the M,-cr re- 
lation would help to constrain theoretical models of bulge 



formation and black hole growth (e.g. ,Silk Rees 1998 



[Burkert SlIk||200T] [Haehnelt Kauffmann|2000[ and oth 



In inactive galaxies, the evidence for the existence of 
black holes and their masses comes from gravitational ef- 
fects on the dynamics of stars inside the black hole's sphere 
of influence. Since the radius of the sphere of influence scales 
with the black hole mass M,, high resolution observations 
are needed to detect SMBHs in the low-mass regime, even 
for the most nearby galaxies. Further difficulties arise from 
the presence of dust in many of these galaxies, particularly 
in discs, which enforces observations in the infrared. Up to 
now the low-mass regime (cr < 120 km s""*^) is sparsely sam- 
pled with only three dynamical black hole mass measure- 
ments (Milky Way, Schodel et al. 2002; M32, Verolme et^aL] 



2002[ N GC 7457|Gebhardt et al.r2003[) and some upp er lim- 



its (e.g. Gebhardt et al. 2001 jValluri et al. 2005). Since 



the near-infrared integral-field spectrograph SINFONI be- 
came operational (Eisenhauer et al. 2003a Bonnet et al. 



2004), it is now possible to detect low-mass black holes in 



dust-obscured galaxies at a spatial resolution close to that 
of HST. 



Based on observations at the European Southern Observatory 
VLT (075.B-0236) 

t Based on observations made with the Advanced Camera 
for Surveys on board the NASA/ESA Hubble Space Telescope 
(GO Proposals 9401), obtained from the ESO/ST-ECF Science 
Archive Facility. 
:|: E-mail: nnowak@mpe.mpg.de 



NGC 4486a is a low-luminosity elliptical galaxy in the 
Virgo cluster, close to M87. It contains an almost edge- 



on nuclear disc of stars and dust (Kormendy et al. 2005). 
The bright star 2.^" away from the centre makes it im- 
possible to obtain undisturbed spectra with conventional 
ground-based longslit spectroscopy. However, it is one of 
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the extremely rare cases, where an inactive, low- luminosity 
galaxy can be observed at diffraction limited resolution us- 
ing adaptive optics with a natural guide star (NGS). This 
feature made NGC 4486a one of the most attractive targets 
during the years between the commissioning of SINFONI 
and similar instruments and the installation of laser guide 
stars (LGS). NGC 4486a is the first of our sample of galax- 
ies observed or planned to be observed using near-infrared 
integral-field spectroscopy with the goal to tighten the slope 
of the M,-cr relation in the low-cr regime (< 120 km s~^) 
and for pseudobulge galaxies. 

This paper is structured as follows: In §2 we present the 
observations and the data reduction. The derivation of the 
kinematics and the photometry is described in §3 and §4 and 
the dynamical modelling procedure is explained in §5. In §6 
the results are presented and discussed and conclusions are 
drawn. 



2 OBSERVATIONS AND DATA REDUCTION 

NGC 4486a was observed on April 5 and 6, 2005, as part 



of guaranteed time observations with SINFONI ( Eisenhauer 
et al.|2003a| [Bonnet et al.|2004| ) at the VLT UT4. SINFONI 



consists of the integral-field spectrograph SPIFFI (spec- 
trometer for infrared faint field imaging; [Eisenhauer et ah] 
|2003bJ and the curvature-sensing adaptive optics (AO) mod- 
ule MACAO dBonnet et al.|2003t . We used the K band grat- 
ing (1.95-2.45 /xm) and the 3'' x 3'' field of view (0.05'' x 0.1'' 
px~^). The bright {R^ 11 mag) star located ^ 2.5" south- 
west of the nucleus was used for the AO correction. The 
seeing indicated by the optical seeing monitor was between 
0.6" and 0.85", resulting in a near-infrared seeing better 
than 0.7", which could be improved by the AO mod- 
ule to reach a resolution of O.l" 25% Strehl). For the 
chosen configuration SINFONI delivers a nominal FWHM 
spectral resolution of ~ 4400. In total 14 on-source and 7 
sky exposures of 600s each were taken in series of "object- 
sky-object" cycles, dithered by up to 0.2". During the last 
observation block the visual seeing suddenly increased above 
1.0". Therefore, to keep the spatial resolution as high as pos- 
sible, only exposures with a seeing < 0.85" were considered 
for the derivation of the kinematics. In order to derive the 
point-spread function (PSF) an exposure of the AO star was 
taken regularly. 



The SINFONI data reduction package spred ( Schreiber 
2006^ was used to reduce the 



et al. 2004 Abuter et al. 



data. It includes all common reduction steps necessary for 
near-infrared data plus routines to reconstruct the three- 
dimensional data cubes. After subtracting the sky frames 
from the object frames, the data were flatfielded, corrected 
for bad pixels, for distortion and wavelength calibrated us- 
ing a Ne/Ar lamp frame. The wavelength calibration was 
corrected using night-sky lines if necessary. Then the three- 
dimensional data cubes were reconstructed and corrected for 
atmospheric absorption using the B9.5V star Hip 059503. As 
a final step all data cubes were averaged together to produce 
the final data cube. Fig. [l] shows the resulting SINFONI im- 
age (collapsed cube) of NGC 4486a. The data of the telluric 
and the PSF stars were reduced likewise. 

As PSF we take the combined and normalised image 
of the PSF star exposures (see Fig.|2|. Its core can be rea- 




Figure 1. SINFONI image of NGC 4486a. 




Figure 2. Left: The two-dimensional PSF. Right: x- and y- 
profiles of the PSF. 



sonably well fitted by a Gaussian profile with a FWHM of 
100 mas (7.7 pc). The diameter of the sphere of influence 
of the black hole can be roughly estimated to 130 mas using 
the M,-cr relation of Tremaine et al. (2002| and is therefore 
resolved. 



3 KINEMATICS 

The kinematic information was extracted using the max- 



imum penalised likelihood (MPL) technique of Gebhardt 



et al. |2000a), which obtains non-parametric line-of-sight ve- 
locity distributions (LOSVDs). As kinematic template stars 
we use six KO to MO stars which were observed during com- 
missioning and our GTO observations in 2005 with SIN- 
FONI using the same configuration as for NGC 4486a. Both 
galaxy and template spectra were continuum-normalised. 
An initial binned velocity profile is convolved with a lin- 
ear combination of the template spectra and the residuals 
of the resulting spectrum to the observed galaxy spectrum 
are calculated. The velocity profile is then changed succes- 
sively and the weights of the templates are adjusted in order 
to optimise the fit to the observed spectrum by minimiz- 
ing the function Xp = + <^^5 where a is the smoothing 
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parameter that determines the level of regularisation, and 
the penalty function V is the integral of the square of the 
second derivative of the LOSVD. We fitted only the first 
two band heads CO2-0 and C03-1. The higher-order band 
heads are strongly disturbed by residual atmospheric fea- 
tures. At wavelength A < 2.29 /im the absorption lines are 
weak and cannot be fitted very well by the templates. 

The uncertainties on the velocity profiles were estimated 
using Monte Carlo simulations (G ebhardt et al. 2QQQaJ . A 
galaxy spectrum is created by convolving the template spec- 
trum with the measured LOSVD. Then 100 realisations of 
that initial galaxy spectrum are created by adding appro- 
priate Gaussian noise. The LOSVDs of each realisation are 
determined and used to specify the confidence intervals. 

In order to test the performance of the method on our 
SINFONI data and to find the best fitting parameters we 
performed Monte Carlo simulations on a large set of model 
galaxy spectra. These were created from stellar template 
spectra by convolving them with both Gaussian and non- 
Gaussian LOSVDs and by adding different amounts of noise. 
We found that the reconstructed LOSVDs resemble the in- 
put LOSVDs very well if the smoothing parameter a is cho- 
sen adequately ( |Merritt|1997| [Joseph et al.|200l| ). The best 
choice of a solely depends on the signal-to- noise ratio (S/N) 
of the data. To maximize the S/N of the data a binning 
scheme with 11 radial and 5 angular bins per quadrant, sim- 
ilar to that used in G ebhardt et al.| ( [2003 ), was chosen. The 
centres of the angular bins are at latitudes 19 = 5.8°, 17.6°, 
30.2°, 45.0° and 71.6° from the major to the minor axis. The 
bins are not overlapping, but spatial resolution elements at 
the border between bins may be divided into parts where 
each part is counted to a different bin. The spectra within 
each bin were averaged with weigths according to their share 
in the bin. The radial binning scheme ensures that an ade- 
quate S/N level comparable to that of the central spectrum 
(S/N?^ 40) is maintained at all radii at the cost of spatial 
resolution outside the central region. 

The resulting two-dimensional kinematics {v, a, /is, h^) 
is presented in Fig. |3] It illustrates the superposition of the 
kinematics of the two distinct components in NGC 4486a - 
the disc and the bulge. Whereas the velocity map shows a 
regular rotation pattern, the cold stellar disc can be clearly 
distinguished from the surrounding hotter bulge in the ve- 
locity dispersion map. The velocity dispersion of the disc is 
~ 20 — 30 km s~^ smaller than that of the bulge. Asym- 
metric and symmetric deviations from a Gaussian velocity 
profile are quantified by the higher-order Gauss-Hermite co- 
efficients /i3 and /i4 ( [Ger hard 1993! [van der Mar el & Franx 



19931 . Fig. [I] shows the kinematic profiles of NGC 4486a 



along the major axis at angles — +5.8° and 6 — —5.8°. 
The +5.8°-profile agrees very well with the adjacent —5.8°- 
profile within the error bars. When comparing the profiles at 
negative radii with the profiles at positive radii slight asym- 
metries can be seen (especially in a), but as the errors are 
relatively large, deviations from axisymmetry are small. 

The star at 2.5^^ from the centre, which heavily dilutes 
optical spectra of NGC 4486a taken without AO, does not 
have any significant effect on the kinematics derived here. 
Quadrant 2 is affected most, as it is located in the direction 
of the star. The fraction of inshining light from the star is 
only about 13% at 0.7^^ from the centre of the galaxy in the 
direction of the star (the outermost point covered by the 




Figure 3. Stellar kinematic fields {v, cr, h^, /i4) of NGC 4486a. 
The binning scheme and the numbering of the quadrants are over- 
laid over the velocity map and the isophotes are overlaid over the 
other maps for comparison. 



exposures of the PSF star) due to the narrow PSF. In addi- 
tion the spectrum of the star shows neither the strong CO 
absorption nor other spectral features in that wavelength 
region. 



4 IMAGING 

To derive the black hole mass in NGC 4486a, it is essential to 
determine the gravitational potential made up by the stellar 
component by deprojecting the surface brightness distribu- 
tion. As NGC 4486a consists both kinematically and photo- 
metrically of two components with possibly different mass- 
to- light ratios T, we deproject bulge and disc separately. 

To decompose the two components, we considered 
the HST images in the broad-band F850LP filter, with 2 
ACS/WFC pointings of 560 seconds exposure each. The two 
dithers have no shift in spatial coordinates. The data were 
reduced by the ST-ECF On-The-Fly Recalibration system, 
see jhttp : //archive . eso . org/ arch ive/hst[ for detailed in- 
formation. 

Moreover, we use the GALFIT package ( Peng et al.|2002 ) 
to fit PSF convolved analytic profiles to the two-dimensional 
surface brightness of the galaxy. The code determines the 
best fit by comparing the convolved models with the science 
data using a Levenberg-Marquardt downhill gradient algo- 
rithm to minimize the of the fit. The saturated star close 
to the galaxy centre has been masked out from the mod- 
elling. The observing strategy, i.e. the adopted no spatial 
shift between the two dithers, has allowed us to obtain a 
careful description of the PSF by using the TiNYTHvQcode. 

We modelled the galaxy light with a double Sersicj 
( 1968| ) model with indices n = 2.19 for the bulge and 



^ Ihttp : //www. stsci . edu/sof tware/tinytim/tinytim 
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Figure 4. Stellar kinematic profiles (-u, cr, /i3, h/C) of NGC 4486a 
along 6 = ±5.8° of the major axis. The corresponding quadrants 
crossed by the profiles are marked for easy identification in Fig. 
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Figure 5. Surface brightness profile for NGC 4486a. At radii left 
of the dashed line we used the decomposed bulge and disc profiles 
for the deprojection, right of the dashed line we used the profile 
of Gebhardt et al. (in prep.) (scaled to ACS z-band). 



Table 1. Results obtained for the four quadrants separately with 
a global mass-to-light ratio (90% C.L.). The numbering of the 
quadrants can be inferred from Fig. [3] 



Quadrant 


M. (10'^ Mo) 


T 


1 


4 0+0.7 


3 4+0.6 

'^•^-0.6 


2 


lQ+1.1 
-^•^-0.2 


o f.+0.2 


3 




o O+0.6 
"^•^-0.1 


4 


-•-•^-1.0 


4 4+0.3 


1 — 4 averaged 




4Q+0.1 



n = 1.67 for the disc. In Fig. [5] we show the surface bright- 
ness profiles of bulge, disc, bulge+disc, and for compari- 
son the SINFONI surface brightness profile and the un- 
decomposed profile of Gebhardt et al. (in prep.), which was 
derived from HST and CFHT imaging in several bands. This 
profile agrees very well with our combined profile and we use 
it, scaled to match our profile, at radii > 11.5^' where the 
bulge strongly dominates. 

Bulge and disc were then deprojected separately using 



the program of Magorrian ( 1999 ) under the assumption that 



both components are edge-on and axisymmetric. The stellar 
mass density then can be modelled as in Davies et al. ( 2006 ) 
via = Tbuige-z^buige + Tdisc-z^disc, whcrc ly is the luminosity 
density obtained from the deprojection and the ACS z-band 
mass-to-light ratio T is assumed to be constant with radius 
for both components. 



5 SCHWARZSCHILD MODELLING 

The mass of the black hole in NGC 4486a was determined 



based on the Schwarzschild (1979) orbit superposition tech- 



nique, using the code of Gebhardt et al. (2000a 2003) in 



the version of Thomas et al. (2004). It comprises the usual 



steps: (1) Calculation of a potential with a trial black hole of 
mass M, and a stellar mass density p*. (2) A representative 
set of orbits is run in this potential and an orbit superposi- 
tion that best matches the observational constraints is con- 
structed. (3) Repetition of the first two steps with different 
values for Tbuige, Tdisc and M, until the eligible parameter 
space is systematically sampled. The best-fitting parameters 
then follow from a x^-analysis. 

The models are calculated on the grid with 11 radial 
and 5 angular bins per quadrant as described above (cf. Fig. 

i- 

Our orbit libraries contain 2 x 7000 orbits. The lumi- 
nosity density is a boundary condition and hence exactly 
reproduced. The 11 x 5 LOSVDs are binned into 17 veloc- 
ity bins each and then fitted directly, not the parametrized 
moments. We limit the parameter space for the values of T 
by considering the population synthesis model of |Maraston| 
( |1998|[2005t , which gives us T < 5 for the z-band. 

Special care was taken when implementing the PSF. 
Due to its special shape with the narrow core and the broad 
wings (cf. Fig. [2]) the PSF was not fitted, rather the two- 
dimensional image of the star was directly used for convolv- 
ing our models. 
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Figure 6. Axq = - X^in ^ function of (top) Tbuige and 
Tdisc5 minimized over M,; (bottom) M, and (^Tbuige + T^^disc) /2. 
The black points are the models we calculated and the coloured 
regions are the (unsmoothed) confidence intervals for two degrees 
of freedom. The best-fitting model is marked with a white circle. 



6 RESULTS 

A big advantage of integral-field data compared to longslit 
data is that we can check the assumption of axisymmetry by 
comparing the kinematics of the four quadrants and quantify 
the effect of possible deviations by modelling each quadrant 
separately. We do not find major differences in the kinemat- 
ics of the four quadrants (cf. Figs, [s] [4|. As it takes a large 
amount of computing time to calculate all models with dif- 
ferent mass-to-light ratios for bulge and disc for each quad- 
rant, we used the same T for bulge and disc for the compar- 
ison of the quadrants. In Table [l] the resulting values for M, 
and T are listed. They show that the four quadrants agree 
reasonably well with each other. The only systematically 
deviant point is the M, determination in the first quadrant 
that however has a large error and therefore is compatible 
within the 90% C.L. with the other three quadrants. There- 
fore we symmetrised the LOSVDs by taking for each bin 
the weight- averaged LOSVDs of the four quadrants and the 
corresponding errors. The results of modelling these aver- 
aged LOSVDs are shown in Fig. [g] where Axo = ~ Xmin 
is plotted as a function of (TTbuige + Tdisc) /2 and M. with 
error contours for two degrees of freedom. Tbuige and Tdisc 
anticorrelate such that their sum is approximately constant, 
as shown in the upper part of Fig. |6] A black hole mass 
of (1.25^0^4) X 10^ Mq (90% C.L.) can be fitted with 
Tdisc ~ 2.8 ... 5.2 and Tbuige ~ 2.8 . . . 4.8. This result agrees 
within 90% confidence limit with the results of all quadrants 
shown in Table^ The best-fitting model, obtained with min- 
imal regularisation, is marked with a white circle and has 
a black hole mass M, = 1.25 x IO^Mq and mass-to- light 
ratios Tdisc = 4.0 and Tbuige = 3.6. The difference in x^ 
to the best-fitting model without black hole is 24.14 which 



corresponds to 4.5cr. The total x^ values for the models are 
around 300. Together with the number of observables (11 
radial bins x 5 angular bins x 17 velocity bins) this gives 
a reduced x^ of ~ 0-3. Note, however, that the number of 



observables is in reality smaller due to the smoothing (|Geb- 
hardt et al.]|2000a| ). 



The dynamical mass-to-light ratios of disc and bulge 
agree with an old and metal-rich stellar population ( Maras- 
ton|1998 2005 ). Tdisc tends to be larger than Tbuige which is 
probably due to the presence of dust in the disc. To estimate 
the effect of the dust on the mass-to-light ratio of the disc 
we are using the model ofjPierini et al. (2004|. For the HST- 
F850LP filter and an assumed typical optical depth r < 0.5 
we obtain an attenuation of Ax < 0.36 mag. This translates 
the best-fitting Tdisc = 4.0 to a significantly smaller dust- 
corrected value of > 2.9. Following the models of Maraston| 



( 1998 2005 ) this is in good agreement with an estimated 



2 Gyr younger disc ( [Kormendy et aL||2005| ). 
The significance of the result is illustrated in Fig. [7| It 
shows the x^ difference between the best-fitting model with- 
out black hole and the best-fitting model with black hole 
(^X^ = YllLi (x?,noBH - X?,bh) over all 17 velocity bins) 
for all LOSVDs of the averaged quadrant. The part out- 
side ^2.5 spheres of influence, where the dynamical effect 
of the black hole is negligible, is displayed in a compressed 
way in order to emphasize the important inner part. For 
75% of all bins the model with black hole produces a fit 
to the LOSVD better than the model without black hole. 
The signature of the black hole is imprinted mainly along 
the major axis, where the largest positive Ax^ are found at 
radii r ^ 0.09'' . . .0.22''. At radii r < 0.09" the differences 
between both models are smaller due to the effect of the 
PSF. For the radii with the largest Ax^ along the major 
axis the LOSVD and the fits with and without black hole 
are shown in the left part of Fig. [7| together with the cor- 
responding Ax? as a function of the line-of-sight velocity. 
The differences between the two fits are relatively small in 
absolute terms. However, the model without black hole has 
more stars on the low- velocity wing at a 1.5cr level, failing 
to match fully the measured slightly higher mean velocity 
of the galaxy. Future observations with a higher spatial res- 
olution should be able to probe this difference more clearly. 

The total stellar mass within 1 sphere of influence, 
where the imprint of the black hole is strongest, is M* = 
9.84x10^ Mq. If the additional mass of M. = 1.25x10'^ M© 
was solely composed of stars, this would increase the mass- 
to- light ratio to Tdisc ~ 9.1 (6.6 if we take into account the 
dust-absorption) , a region which is excluded by stellar pop- 
ulation models or at least requires unrealistic high stellar 



The models with black hole become tangentially 
anisotropic in the centre, while the models without black 
hole are close to isotropic. 

Our result is in good agreement with the prediction of 
the M.-cr relation (1 .26^0 35 x lO'^ M© using the result of 



Tremaine et al.|2002 ) and strengthens it in the low-cr regime 
(< 120 km s~ ), where, besides several upper limits, up to 
now only three black hole masses were measured with stellar 
kinematics (Milky Way, [Schodel et al.||2002| M32 



Verolme 



eraL][2002l NGC 7457, [Gebhardt et al.|2003 ). 

NGC 4486a is only the first object in our sample of 
low-mass galaxies under investigation. The laser guide star. 
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Figure 7. (a) difference between the best-fitting model without black hole and the best-fitting model with black hole (Ax^ = 
y^ ^. Ax^ = X^i-i (/^inoBH ~ ^?bh) -'-'^ velocity bins) for all LOSVDs of the averaged quadrant. Bins where the model with 

black hole fits the LOSVD better are plotted in green, the others in orange. The part outside ^2.5 spheres of infiuence is plotted with a 
different scale than the inner part, since in the outer region the dynamical effect of the black hole (and therefore the difference between 
the two models) is negligible. (b,c) For the radii with the largest positive x^ difference in (a) the LOSVD (open circles with error bars, 
normalised as in |Gebhardt et al.|2000a| > and both fits (with black hole, full green line, and without black hole, dashed orange line) are 
shown with the corresponding Ax^ plotted below (d,e). 



which is presently being installed at the VLT UT4, makes 
observations of a large number of appropriate galaxies now 
possible. Therefore we plan to further explore this region 
of the M,-cr relation by observing more galaxies with ve- 
locity dispersions between the resolution limit of SINFONI 
(- 30 km s"^) and - 120 km s"\ 
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